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Abstract 

Within a scaling ansatz formalism plus Thomas-Fermi approximation, we 
investigate the collective excitations of a harmonically trapped boson-fermion 
mixture in the collisionless and hydrodynamic limit at low temperature. Both 
the monopole and quadrupole modes are considered in the presence of spher- 
ical as well as cylindrically symmetric traps. In the spherical traps, the fre- 
quency of monopole mode coincides in the collisionless and hydrodynamic 
regime, suggesting that it might be undamped in all collisional regimes. In 
contrast, for the quadrupole mode, the frequency differs largely in these two 
limits. In particular, we find that in the hydrodynamic regime the quadrupole 
oscillations with equal bosonic and fermionic amplitudes generate an exact 
eigenstate of the system, regardless of the boson-fermion interaction. This 
resembles the Kohn mode for the dipole excitation. We discuss in some detail 
the behavior of monopole and quadrupole modes as a function of boson- 
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fermion coupling at different boson-boson interaction strength. Analytic so- 
lutions valid at weak and medium fermion-boson coupling are also derived 

and discussed. 

PACS numbers:03.75.Fi, 05.30.Fk, 67.60.-g 
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I. INTRODUCTION 



Shortly after the achievement of Bose-Einstein condensation of dilute, magnetically 
trapped alkali atoms [0,0, the investigation of collective excitations in these systems has 
become a very active research field (see for a recent theoretical review). The high accu- 
racy of frequency measurements and the sensitivity of collective phenomena to interaction 
effects makes them good candidates to unravel the dynamical correlation of the many-body 
system. So far, experimental results have been obtained for low-lying collective modes of 
a trapped condensate in a wide temperature regime, including breathing modes [Q, sur- 
face modes p, and the scissors mode p. These experiments have in turn stimulated a 
considerable amount of theoretical work. 

Recently the quantum degenerate regime has also been reached in a magnetically trapped 
Fermi gas 0, and in a mixture of Bose and Fermi particles The latter system is in 

particular interesting since it serves as one typical example in which the intermingled par- 
ticles obey different statistics. Up to date the static property , the phase dia- 
gram and phase separation [|14|,|15|,|16[ , stability conditions ||T7|JT8|| and collective excitations 
2^p5| , |26[| of trapped boson-fermion mixtures have been theoretically investi- 



gated. In a recent experiment, the collapse of a degenerated Fermi gas caused by the strong 
attractive interaction with a Bose-Einstein condensate has been observed in an atomic mix- 
ture of ^°K— ^^Rb and measurements of collective excitations might be available soon 
also in such systems. 

The purpose of the present paper is to study the collective excitations of magnetically 
trapped boson-fermion binary mixtures in two different regimes: a collisionless regime where 
the collision rate is small compared with the frequencies of particle motion in the trap and 
a hydrodynamic (collisional) regime in which collisions are sufficiently strong to ensure local 
thermodynamic equilibrium. From the experimental point of view the temperature T of all 
the realized boson-fermion mixtures is around the Fermi temperature Tp (more precisely, 
T ^ 0.2Tp M), and the systems are possibly in or close to the hydrodynamic regime, since 
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the sympathetic coohng technique used in the experiments usually requires a large boson- 
fermion interaction strength so that the frequent collisions between fermions and bosons 
can foster the local thermal equilibrium and ensure efficient thermalization of the fermionic 
component to reach the quantum degeneracy. Of course, as far as the strongly degenerate 
regime (where T -C Tp) is concerned, the collisions are rare because of Fermi statistic, and 
the systems will be finally in the collisionless regime. 

Several theoretical analysis have already been presented for the collective excitations of a 
spherically trapped boson-fermion mixture. The collisionless modes are considered by a sum- 



rule approach |24] or in the random-phase approximation |25,26|. The collisional collective 



oscillations are discussed by Minguzzi and Tosi [23|, however, limited to the surface modes at 



weak fermion-boson coupling. On the other hand, the homogeneous boson-fermion mixtures 



have also been analytically studied |T9lp0lpT,22|. The repulsion between the Bogoliubov 



phonon mode and zero- sound mode |jT9|] (or Anderson mode 0]), is predicted when the 



degenerated Fermi gas is in the normal collisionless limit (or in the superfluid phase). 

In this paper, we shall analyze systematically the collisionless and hydrodynamic 
monopole and quadrupole modes in the presence of spherical as well as cylindrically sym- 
metric traps [^. We discuss in some detail the behavior of those modes against boson- 
fermion coupling at different boson-boson interaction strength. In the spherical traps, we 
find that the monopole frequency coincides in the collisionless regime and in the hydrody- 
namic one, suggesting that it might be undamped in all collisional regimes. In contrast, for 
the quadrupole mode the frequency differs dramatically in these two limits. In particular, 
in the hydrodynamic regime the quadrupole oscillations with equal bosonic and fermionic 
amplitudes are found to generate an exact eigenstate of the system, resembling the Kohn 
mode for the dipole excitation. Analytic solutions valid at weak and medium fermion-boson 
coupling are also deduced and discussed. 

The content of the paper is as follows. In the next section we derive equations of the low- 
energy collective excitation of a boson-fermion mixture in the Thomas- Fermi approximation 
by means of a scaling ansatz. In Sec. Ill we first briefiy describe the parameters and 



the numerical procedure employed in the present calculation. We then turn to detailed 
discussion of the results obtained by analyzing the dependence of the mode frequencies on 
the boson-fermion coupling and boson-boson interaction strength, the ground-state density 
distributions, and the mixing between bosonic and fermionic collective oscillations. The last 
section is devoted to summary and conclusions. 



II. FORMULATION 

We consider a dilute spin-polarized boson-fermion mixture trapped in a cylindrically 
symmetric harmonic oscillator potential at low temperature. In the semi-classical Thomas- 
Fermi approximation, the normal Fermi gas evolves according to the Boltzmann-Vlasov 
kinetic equation |2^j3^ (see the Eq. (0) below). In the coUisionless regime, the collisions 
are rare and we can safely neglect the collision term (Icoii) that accounts for the damping of 
collective modes. On the opposite of the hydrodynamic regime, the collision term dominates 
and we resort to the Euler equation of motion |31[], which can be deduced from Boltzmann- 



Vlasov kinetic equation under the assumption of local equilibrium for fermions |32|. For the 



bosonic part, we shall apply the same Stringari's hydrodynamic formulation [p3| in both 



regimes, since the dynamics of the condensate is less affected by the collisions 13^ . 



A. coUisionless regime 

According to the Stringari's hydrodynamic description PB[, the low-energy dynamics of 



the trapped bosonic atoms is determined by the equations 

drib 



dt 



+ V (vbUb) = 0, 



"^6-^ + V (^V^„ + ^f,bnb + ^f,jn/ + -mhVftj =0, (1) 

where I^m (r) = \mh {oj\_hP^ + '^26^^) is the cylindrical symmetric confining potential, and 
nf, (r, t), Tij (r, ^) and Vf, (r, t) are the boson, fermion density and velocity field, respectively. 
The mean-field term ghfn^ (r, t) is included to take into account the effect of boson-fermion 



interaction [|T0|. The boson-boson and boson-fermion interaction strength for the pseu- 
dopotentials, Qbt and Qbf, are related to the s-wave scattering lengths a^fe and a^f through 
gbb = 'i'^'h'^abb/mb, gbj = 2Txfi^abf /mbf where mbf = mbUif/ {nib + mj) is the reduced boson- 
fermion mass. In Eq. (|1|), we have already neglected the quantum kinetic energy pressure 
V^A/?Tb in the spirit of the Thomas- Fermi approximation [p[. 



term tt-^ , — 

For the fermionic part in the coUisionless regime, in order to take into account the 
effects of the boson-fermion interactions, we consider the mean-field description based on 
the Boltzmann-Vlasov kinetic equation p9i|30| without the collisional term (Icoii), 



dl 
dt 



df 1 dVl df gbfdrib df 



0, 



(2) 



dr uif dr dvf mj dr dvj 
where /(r, v/,t) is the single particle phase space distribution function for fermions, 
n/(r, t) = / (i^vj/ (r, vj, t) and V/„ (r) = ^mj {uj^^p^ + ujI^z^'^ is the confining potential. 
The last term in the left hand side of Eq. (0) is a Hartree-Fock mean-field term, also 
known as the Vlasov contribution. The fermion-fermion interaction has been neglected as 
the polarized system is considered [BB . 

Without the boson-fermion interaction [gbj = 0), both the Eqs. ([l|) and admit the 
simple scaling solution, i.e., 

1 n / ri \ 



rib (r,t) 



Vbi (r, t) 



1 dbi{t) 



hi{t) dt 



-r,;, 



(3) 



for bosons 31.36 and 



/(r,v/,t) = /o 



Mi) 



,v/(r,t) 



v/i(r,t) =^i(t)vfi —n 



(4) 



for fermions p9| , |30| , |37| . Here and /o are the equilibrium distributions. The dependence 
on time t is entirely contained the six scaling parameters, hi (t) and 7i(t), where i = x, ?/,and 
z. By substituting this solution into Eqs. (|I]) and (^, it is easily to show that, the scaling 



parameters obey the coupled differential equations p6| , p7 | 



0, 



(5) 



kit) Ujbjit) 




7,(t)+u;5(t)7,(t) 



0. 



(6) 



Solutions of Eqs. (|^) and (^) determine the evolution of both the boson and fermion density. 
In particular, the eigenfrequencies of small oscillations with uJi{t) = uJi{0) are the resonance 
frequencies of the collective density (shape) oscillations under the weak perturbation of the 
external field. After linearizing around the equilibrium values bi = 1 and 7^ = 1, one finds, 
in the case of the spherical harmonic traps, the result: UMb = V^^^b, ^Qb = V^^^b and 
UMf = '^^f, ^Qf = Scu/ for the frequencies of the monopole and quadrupole oscillations, 
respectively, which is already well known in the literatures. 

In the presence of the boson-fermion interaction {gi,f 7^ 0), however, the simple scaling 
solution is no longer satisfied at every position r after the substitution. A useful approx- 
imation, in the first order of ghf-, is to assume the scaling form of the solution as a priori, 
and fulfill it on average by integrating over the spatial coordinates. The same strategy has 
been recently used by Guery-Odelin p9[ to investigate the effect of the interaction on the 
collective oscillation of a classical gas in the coUisionless regime and by Menotti et al. to 
study the expansion of an interacting Fermi gas ^0[. In the latter, the authors showed that 
the frequencies of the monopole and quadrupole modes for isotropic traps deduced in this 
approximation coincide with the result derived earlier by using the sum-rule approach. Of 
course, the assumption of the scaling ansatz is only meaningful for the small value of \ghf\- 
It will apparently break down for a large and positive ghf-, at which the phase separation 
occurs. 

In the approximation specified above, we substitute the scaling ansatz Eq. (j^) into 
Stringari's hydrodynamic equations. By setting Ri = ri/bi(t), one finds. 



bi{t)Ri + ul{t)bi{t)Ri + 




(7) 



mbbi{t)Ujlj (t) dRi 
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which reduces to the following form in equilibrium state, 

+ a^Hm , ^MM . 0. (8) 

rrih oRi mh oRi 
The equations for the scaling parameters 6j(t) can be obtained by multiplying Eq. (|^) by 
RiU^ (Ri) on both sides and integrating over the spatial coordinates. Making use of the 
equilibrium properties of the density distribution (||), after some straightforward algebra 
one finds, 



where {Ri)f, = ^ I d^^nl(TV)Rf is the average size of bosons along the i-axis. The last two 
terms in Eq. (|^), linear in Qbf, account for the effects of boson-fermion interaction. 

Analogous procedure can also be applied for the fermionic part The equations for 
the scaling parameters 7j(t) are finally found to take the form. 



9bf 1 frfSK^o^oA. 



'mjNf{R^)f^,UjljJ OR, ' ^S, 



where (fl.^)^ = J. / #Rn»(R)JJ?. 

The coupled set of differential equations @ and ([T0|) is a generalization of Eqs. (j^) and 
(^) in the presence of the boson-fermion coupling. It determines the dynamics of boson- 
fermion mixtures in the collisionless regime as far as the assumption of the simple scaling 
solution is valid. We shall only be interested in the small oscillation around the equilibrium 
state (6i,7i ~ 1) and apply it to study the behavior of monopole and quadrupole modes 
against boson-fermion coupling. In this case, one can simplify the set of differential equations 
by expanding 
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to obtain 



+ E.#^f?^-l)-0, (12) 



and 



hit) Ujbjit) 

'k 



+i:»^^A;^(?-i1 = o. (13) 



where 



R /,3p ^^^(R) pp g^g(R) 



mi,u;l{0)Nt{R^)J dR, ' " dRk ' 

mfUjtj{0)Nf (Rf) J J dRi ^ 



mfu;ff{0)Nj{Rf)^J dR, ' ' dR^ ' 

are the dimensionless parameters proportional to gbf. In the case of a cyhndrical trap, those 
parameters can be reduced to Ba/3, Fa, and Da/3 {c(,P = P or z) in terms of the cyhndrical 
coordinates, whose expressions are given in Appendix. 

By linearizing Eqs. (|l^) and ([T3| ) around fej, 7^ = 1, for each component in p or 2; coordi- 
nate one gets a separated equation and thus obtains four coupled equations. The dispersion 
relation for the frequency of the monopole and quadrupole modes can be determined by the 
condition for existence of nontrivial solutions, that is. 



det 

where the matrix 



0, (14) 
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A, 





For each frequency of modes, the corresponding nontrivial solution, denoted by {Abi,Ab2,Afi,Af2), 
gives the amphtude of the small density oscillations. As the system is composed of two kinds 
of particle, one would expect an emergence of two type of collective oscillations for each mul- 
tipole. We thus define the mixing angle 



to characterize the degree of mixing between bosonic and fermionic collective motions. As 
a limiting case, 9 = 11/2 ot corresponds to the purely (decoupled) bosonic or fermionic 
oscillations, respectively. 

In the special case of spherical traps, the monopole and quadrupole modes are decoupled. 
As shown in the Appendix, the dimensionless parameters Bajs {Fa, Dap) can further be 
reduced to a single value B (F, D), 



Accordingly, the dispersion relation (|T4D for the frequency of modes takes the following 
simple form. 




(15) 




(16) 




(18) 



(17) 
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where the suffix M, Q denote monopole and quadrupole modes, respectively. 



B. hydrodynamic regime 

In the hydrodynamic regime, one assumes that local equilibrium has been established 
for the ground-state density proffies and it is maintained during dynamic fluctuations of the 
particle densities. A useful description of the dynamics of the Fermi gas in this regime is 



based on the Euler equation of motion |31,3 



dt 

dvfi 1 dP(r,t) ^ dvfi d \ 

where the pressure P{r,t) = | — — ^/''^(i')^) the local density approximation . 
Without the boson-fermion interaction {gbf = 0), Eq. (|T^) still admits a simple scaling 
], i.e., 

1 



V (YfUf) = 0, 



solution [31 



n/(r, t) =— —n 



nn,{t) '\-ii(t)J 

which leads the following equations for the scaling parameters 7i(t): 

m + u^UtMt) ,2/3 = 0- (21) 

In the presence of a nonzero g^j, one may still expect the validity of such scaling ansatz in 
the weakly coupled limit, although there is no verification by the other method (Note that 
the sum- rule approach is only applicable in the collisionless regime). Along the same line 
as in the previous subsection, after substituting the scaling ansatz (^OD into equation (p!9D 
and taking the moment with respect to rf, one ultimately finds the differential equations 
satisfied by 7j(t), 
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9bf 



It 



9bf 



li 



V3./ rf^R^i?.-;(R) = o 



(22) 



which differ from Eqs. ([lOD for the coUisionless regime. 

Eqs. (H) and (p2|) can be linearized around the equilibrium state and combined to yield 
the determinant, 



det 



0, 



(23) 



A, 



(4 




(1 


- Bpz)u:l, 




Bpp^\b 






(2 


- Bzp)u'^f, 


(3 


- B^^)u^f^ 




Bzpi^zb 








DppUj'if 




Dpzi^± f 


{- 

V 3 


3 ^Pp)^±f 


{- - 

V3 


-y - Dpz)uj'j_j 










{- 


2-^^ n V ,2 


- 

V3 





which gives rise to the dispersion relation in the hydrodynamic regime. 

Finally, the frequency of the monopole and quadrupole modes for the spherical traps is 
obtained by rewriting Ba/3 {Fa, D^p) in terms of B (F, D), 

o 1 



M 



2 

± 



{[(5-5; 



{A-F- D)ujj 



(5 - B)iut - (4 - F - D)iuj + ABDiutujj 



1/2 ~ 



u;l = -{[{5-B)u;', + {5~D)u} 



5 
± 



(5 - B)iul - (5 - D)iu'j 



+ ABDujIuj) 



1/2^ 



(24) 



(25) 



Eqs. (p!7| ) and ( p^ explicitly show that the frequency of monopole oscillations coincides in the 
coUisionless and hydrodynamic regime. This fact is a reminiscent of properties of a classical 
gas confined in an isotropic traps {uo^ = ujy = uo^ = oJq) [^,0, in which the monopole 



oscillation of frequency uj = 2ujo is an exact undamped solution of the full Boltzmann 
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equation. Analogously, in the boson-fermion mixtures the monopole excitation might also 
be undamped in all coUisional regimes from the coUisionless to the hydrodynamic one. 

At the end of this section, we briefly mention the whole process of the numerical cal- 
culations that consists of three stages. First, one has to find the equilibrium ground-state 
densities at low temperature, which approximately satisfy the following coupled equations 
in the Thomas- Fermi approximation |[10[| , 



Vho (P> + dbbrib (p, z) + gbfU^f (p, z) = /it, 
(evr^nj (p, z)f + Vl (p, z) + g.^nl (p, z) = fif, (26) 

where pbj is the chemical potential. Then one computes the dimensionless parameters -B^/j, 
Fa, and Dap {a,l3 = p or z), and finally, one solves the Eqs. (|14D and (^31) (in the case of 
spherical traps, one explicitly uses Eqs. dT^); ( P^ and (^)) to obtain the frequency of the 
monopole and quadrupole modes. The mixing angle for each mode is also simultaneously 
calculated. 

III. RESULT 

In this work, we have performed a numerical calculation for Nh = Nf = 10®, where the 
number of bosons and fermions is large enough to ensure the validity of the Thomas- Fermi 
approximation, i.e., Ni,aiih/a\ ^ 1 and Nf ^ 1 [|l^,|l^. We take the harmonic oscillator 



length a\ = w — - — and huj±h as units, and define the scaled dimensionless variables: the 
coordinates p = p/a\, z = z/a^, boson/fermion densities n^j = n^j {o.±^ , interactions 
strength g^i, = gb^/ huj±b {o.±) — 4vrafeft/a5_, and chemical potentials jlbj = Pbj/ {h^ib}- 
We also introduce the quantities a = mf/rrib, [3 = u!±f/uj±b, A = u!zb/^±b = ^zf/^^f, 
and K = gbf/gbb to parameterize the different mass of the two components, anisotropy of 
traps, and boson-fermion coupling relative to the boson-boson interaction. The constraint 
= 1 is always satisfied since in experiments both bosons and fermions experience the 
same trapping potential. In the scaled units, the coupled Thomas- Fermi equations ( PB] ) take 
the form, 
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It is convenient to obtain the solutions to Eq. ( P7D by iterative insertion of one density 
distribution in the other equation and numerically searching for the chemical potential jib 
and fif yielding the desired number of particles. 

We shall investigate the behavior of monopole and quadrupole modes against boson- 
fermion coupling n for three typical values of Qbt ■ First of all, we consider the relevant 
parameters for a boson- fermion mixture composed of (fermion) and ^^Rb (boson), which 
has been recent realized by the LENS group In order to emphasize the interplay of 
collective modes of bosons and fermions due to the nonzero gbf (the degree of mixing will 
be maximum if theirs bare mode frequencies are close to each other), we will consider the 
same mass (m^ = rrif = m) and trapping frequency {uJb = ujf = ujq) for bosons and fermions 
in most cases, although the realistic mass of ^'^Rb is about two times larger than that of 
^°K. As in the experiment, we take the radial harmonic frequency of uj±b = 27r x 216 s^^ 
for ^'^Rb and the boson-boson s-wave scattering length of abb = llOao = 5.9 nm, which 
gives the rescaled interaction strength cjbb = 0.1. In the ground state the fermions have a 
much broader distribution than bosons because of the Pauli principle and the bosons are 
completely immersed in the Fermi sea [^]. Secondly, as an opposite limit, we consider the 
case in which the fermions and boson have approximately the same radius and significantly 
overlap with each other. Within the Thomas-Fermi approximation, at k = 0, the radius 
of the Bose condensate and zero-temperature Fermi gas in the scaled units are given by 
Tb = (IbNbgbb/'i'^Y^^ and rj = (ASNfY^^, respectively. Equating these two numbers we get 
the constraint: Qbb = 2.11. Finally, we take c/bb = 0.5 for the intermediate regime. It should 
be note that in experiments the interaction strength cjbb can be controlled by using Feshbach 

m. 



resonances 



In the section IIIA, we briefly estimate the criterion for establishing the hydrodynamic 
regime. In the next sections IIIB and IIIC, we analyze the collective modes in the collisionless 
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and hydrodynamic regime for a spherical trap. The results for a cylindrically symmetric trap 
will be presented in section HID. 



A. the criterion to establish the hydrodynamic regime 

A hydrodynamic regime is established in the low temperature alkali vapor when the 
inequality 

a;T < 1 (28) 

holds, where r being the collision time for incoherent scattering of fermions against the 
condensate and uj being on the scale of the trap frequency (a; — uJf) for the low-lying modes. 
At low temperature, the dominate collision procedure comes from the scattering between a 
fermion and a condensate boson, which generates another fermion and a Bogoliubov quasi- 
particle. In this procedure, the mean velocity of the condensate boson is negligible relative 
to that of the fermions. A naive estimate for the collision time can thus be written as 

T-^^n,(AT:alf)vF(^^^ , (29) 

where {Anal^ vp is the classical coUisional frequency and the factor (^) results from 



the Pauli blocking. By setting = iV(,/(47rr^/3) and taking vp = \j2iip/mf for a spherical 
trap, we approximately have 

_L.3x2-.6V.x(J|-)"^-^ra^,(^)\ (30, 

where o;,/ = abf/aho is the s-wave boson-fermion scattering length in the scaled units. 

For illustrative purposes we again consider the boson-fermion mixture of and ^^Rb 
studied by the LENS group. Prom the known values of the ^°K-^'''Rb scattering length 
ttbf — SOOoo, and = 0.1, we have 

1 / T \^ 

29 X — , (31) 



OUfT \Tp 
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for A*";, = Nf = 10^. As anticipated earlier in the introduction, the temperature of order 
T ~ 0.5Tp would suffice to verify the inequality ( pHD with uj ~ Uf, and therefore to reach 
the hydrodynamic regime. 

In the next subsections, we shall consider the behavior of the collective modes against 
the boson-fermion interaction strength, rather than the temperature. It should be reminded 
that for a fixed temperature, our results for the hydrodynamic regime (or the collisionless 
regime ) are only valid at \gbf \ ^ g^j (or \ghf \ ^ g^jr), where g^jr can be roughly determined 
from Eq. (ID. 



B. collisionless modes in a spherical trap 

Figures 1 and 2, respectively, show the frequencies and mixing angles of the monopole 
and quadrupole modes as a function of n for the case in which the bosons and fermions 
have the same mass (m^ = rrif = m) and trapping frequency (tUf, = ujf = ujq). As we 
have anticipated in the last section, there are two types of collective oscillations for each 
multipole. For clarity, we plot the lower and higher frequency modes by the solid and dashed 
lines, respectively. 

The collisionless collective modes for the weak boson-boson interaction has been inves- 
tigated earlier by a sum-rule approach in Ref. As shown in the figures la and 2a, the 
result for gtb = 0.1 agrees well with that obtained by the sum-rule approach in the whole 
regime of n if we use the same parameters (see, for example, the figures 2a and 2b in Ref. 



2J]). This excellent agreement in some sense justifies our assumption of the scaling form of 
solutions to Eqs. (|l]) and (0). We believe the same is true in the hydrodynamic regime. 

The most remarkable feature in figures 1 and 2 is the existence of a specific value Kc ^ 
, at which the collective oscillation of each mode becomes purely bosonic or fermionic. Kc 
coincides with the critical values of phase separation for the strong boson-boson interaction 
and becomes unity for the weak or medium boson-boson interaction. The existence of Kc 
can be readily understood from the distribution of boson and fermion densities. As shown 
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by M0lmer [^], for a positive boson-fermion interaction, the fermions are squeezed out the 



center. As k increases, they will eventually form a shell-like distribution around the surface 
of bosons for K ^ 1 and will be completely pushed away from the center at a critical value 
where the phase separation occurs. For the weak or medium boson-boson interaction, as 
shown in the figure 3a, precisely at k = 1 the fermions experience a constant potential in 
the region occupied by bosons and therefore uniformly distributed there. As a result, the 
parameters B and D defined in Eqs. (|16|) will vanish and consequently the bosonic and 
fermion part in the determinant (|T^) will be completely decoupled. Therefore the collective 
oscillation of modes will be purely bosonic or fermionic. Analogous mechanics works for the 
case of strong boson-boson interaction, where B and D will be zero at the critical value of 
phase separation that is smaller than unity. Note, however, that in this case our assumption 
of the scaling ansatz will apparently break down once k, > and the phase separation will 
lead to a decoupling of the collective modes. 

An immediate application of the above observation is that we can derive an analytic 
expression for the frequency of each mode in the weak or moderately strong boson-boson 
interaction in which = 1. As illustrated in figure 3b, for small values of k, the B, F, 
and D can be well approximated by bQK,{l — k), /qk, and dQK,{l — k), respectively, where 
bo = {dB / dK) , /o = {dF / dK) , and d^ = {dD / dn) ^^q. The form of B and D follows the 
fact that they have to vanish at both k = and k = 1. In the Thomas- Fermi approximation, 
one may obtain, 

224 fNf\ , /• , / , ,Ni/2 





fo = -^-^x'Jy'{l-xVf"dy, 



1 



do = +^x'^ Jy'{l-xVy^'dy, (32) 





where x = ^ = (^)^^^ (15A/'fe^fe;,/47r)^/V(48iV/)i/6 ^ i is the ratio of the radius of the 
Bose condensate and zero-temperature Fermi gas. The frequencies obtained by combining 
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Eqs. ([TTD, ([IHD and (0) are plotted in figures 1 and 2 by thin lines. We find that it is in a 
good agreement with the full numerical calculations for a wide regime of k. 

Below we discuss the behavior of the frequencies Ua of each mode by defining three 
regions of k: (I) k < 0, (II) ^ k < Kc, and (III) k,c ^ 

(a) Monopole. For a noninteracting boson-fermion system the low-lying monopole mode 
is the fermionic oscillation with frequency uj^j = 2ujq, while the higher mode is the bosonic 
one with uj^j = V^luq in the Thomas- Fermi approximation. Around k = 0, one may obtain 
cuf'j ^ 2ujq (l — '^°g-^° ft:^ and u^j ~ V^ujq (l — fo'^)- '^^^ essential features of the monopole 
mode as a function of n can be summarized as follows, (z) The curve for the high-lying 
bosonic mode seems like a parabola. Indeed, the frequency for the bosonic mode is always 
found to be degenerate at k = and k, = Kc- In region II, the frequency varies slowly 
against n and is slightly smaller than the value of V^ujq in the noninteracting limit. Here, the 
boson density distribution expands slightly compared with that of the noninteracting boson- 
fermion system due to the weakly repulsive boson-fermion interaction and simultaneously 
the bosons experience a weaker effective confinement. In the region I and III, the situation is 
quite different: the bosons are heavily compressed by either the attractive or strong repulsive 
boson-fermion interaction. This strong confinement leads to a steep rise of frequency with 
increasing \k\. (ii) The behavior of the frequency for low-lying fermionic mode against k 
is more complicated. As the boson-boson interaction strength gi,h increases, the sign of the 
derivative of the curve at k = changes from positive to negative. At g^b = 2.11 a. large dip 
appears in the region II. On the other hand, at a large negative values of k, as pointed out 



by Miyakawa et al. [^, we find a sharp decrease of the frequency towards the instability 
point of the ground state. {Hi) Finally, at large values of |k| the mixing angle for both 
the low-lying and high-lying modes becomes close to 7r/4, suggesting that the bosons and 
fermions are highly correlated in the collective oscillations. The degree of mixing is enhanced 
as one increases gbb- 

(h) Quadrupole. For the quadrupole excitation (figure 2), the lower and higher energy 
mode becomes bosonic and fermionic oscillation, respectively. To the first order of k the 



frequencies of the lower and the higher quadrupole modes are given by Uq ~ V^ujq (l — j^hzj 
and lJq ~ 2ujq (l + ^'^°2o'^° ^) • For those modes, similar mechanisms as for the monopole mode 
are still at work concerning the dependence on k. However, the role of the boson-fermion 
interaction is much reduced compared with the monopole case as seen by the factor 1/5 in 
Eq. (jT^), which reflects that the quadrupole oscillation has 5 different components [^. In 
addition, the behavior of the high-lying fermionic mode is somewhat simpler. The frequency 
is always decreases around k = as /t increases. 



C. hydrodynamic modes in a spherical trap 

In section II, we have analytically showed that in a spherical trap the frequency of 
monopole oscillations coincides in the coUisionless and hydrodynamic regime. Here we find 
a dramatic difference for the quadrupole mode. In the figure 4, we plot the frequencies 
and mixing angles of the quadrupole modes in the hydrodynamic regime against k. The 
parameters are the same as that in figures 1 and 2. Compared with the result for the 
quadrupole mode in the coUisionless regime (figure 2), an interesting feature emerges: The 
frequency of the low-lying mode (high-lying mode) is always fixed to V^coq in the region 
I and III (region II), independently of the value of k, and the corresponding mixing angle 
is exactly 7r/4 even around k = 0. This strongly suggests that in this case the collective 
oscillation with equal bosonic and fermionic amplitudes generates an exact eigenstate of the 
system, regardless of the boson-fermion interaction. It resembles the Kohn mode in the 
isotropic harmonic traps. Indeed, the behavior of the frequencies shown in figure 4 is quite 
similar to that of a dipole mode in the coUisionless limit (see, for example, the figure 2c in 
Ref p3). The above feature can also be explained explicitly from Eq. (pS]). For the case 



considered here, Ub = ujf = ooq, "we have ujq = y/2uJo {1 — [{B + D) ± \B + D\] /lO}^''^, and 
thus one branch of ujq will always be v^Wq. 

The strong mixing of the bosonic and fermionic oscillation stated above in fact arises from 
the "on resonance" condition, that is, the frequency of the bosonic and fermionic quadrupole 
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modes in hydrodynamic regime is degenerate in the noninteracting limit of k — 0. As shown 
in figure 5, once one moves away from the "on resonance" condition by changing mf/mb, 
the degree of mixing will be much reduced. 

D. a cylindrically symmetric trap 

In this subsection, we consider a cylindrical trap that is more relevant to the experiment. 
Figures 6 and 7, respectively, display the frequencies of each mode for a cigar-shaped (A = 
0.5) and disk-shaped (A = 2.0) trap in the coUisionless (thick lines) and hydrodynamic 
(thin lines) regime. As the monopole and quadrupole modes are coupled away from A = 1, 
we denote the two higher and two lower modes as quasi-monopole and quasi- quadrupole 
ones, respectively. In the literature, the former one is also called as transverse breathing 
mode in the limit of A ^ 0. As one may expect, in the cylindrical trap the frequency of 
quasi-monopole mode in the coUisionless and hydrodynamic regime is no longer degenerate. 

(a) a cigar-shaped trap. For the quasi-monopole excitation (figure 6a), the lower and 
higher energy mode in both regimes are fermionic and bosonic oscillations. In region I and 
III, the frequencies for the high-lying mode in the coUisionless and hydrodynamic regime 
are almost the same and they only differ sightly in region II. In contrast, for the low-lying 
mode the frequencies in two regimes have significant differences. In particular, as k decreases 
towards the instability point of the ground state, the frequency in the coUisionless regime 
rises up steeply, while the one in the hydrodynamic regime shows a sharp decrease. For the 
quasi-quadrupole excitation (figure 6b), on the other hand, the frequencies for each mode 
bears a lot of similarity as that in the spherical case. We thus don't discuss them further. 

(b) a disk-shaped trap. In this case (figure 7), the lower and higher energy mode for both 
quasi-monopole and quasi-quadrupole excitations in the coUisionless regime are bosonic and 
fermionic oscillations, respectively. The opposite is true in the hydrodynamic regime. Most 
interestingly, the soften of the mode frequency towards the instability point of the ground 
state now appears in the low-lying mode for the quasi-quadrupole excitations, rather than 
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the monopole one. 



IV. SUMMARY AND DISCUSSION 

Ultracold boson-fermion mixtures of alkali atoms have recently been the subject of in- 
tensive experimental research. As an important tool to characterize the behaviour of this 
kind of many-body system, the investigation of collective oscillations will be of particular 
interest. In this paper, with the help of Thomas- Fermi approximation and a scaling so- 
lution we have studied the behavior of monopole and quadrupole excitations against the 
boson-fermion interaction in two limiting cases: the collisionless and hydrodynamic regime. 
For a spherical trap, the frequency of monopole mode is identical in both regimes, analo- 
gous to that of a classical gas, which is undamped in all collisional regimes. In contrast, 
the frequency of quadrupole mode differs largely in these two limits. Most interestingly, in 
the case of same trapping frequency for the two components {oob = ^f) and hydrodynamic 
regime, the quadrupole oscillations with equal bosonic and fermionic amplitudes are found 
to generate an exact eigenstate of the system, regardless of the boson-fermion interaction. 
It indeed resembles the Kohn mode for the dipole excitation. 

While we have restricted the discussion to the collisionless and hydrodynamic regime, it 
should be explicitly remarked from the experimental point of view that it is more interest- 
ing to investigate the crossover between these two limits , which might be realized in 
the experiment by changing the temperature or controlling the particle-particle interaction 
through Feshbach resonances. In theory, such crossover might be investigated by adding a 
collision term {Icou) to Eq. (||). 

We are aware of that the above results are based on the assumption of the simple scaling 
solution. Its validity has been partly justified by the good agreement between our results and 
that obtained by a sum-rule approach. However, in case of large boson-fermion interaction. 



the spectrum of lowest collective excitations might be fragmented p5| , p6| , |i3[] . As a result. 



both the sum-rule approach and the approximation of simple scaling solution will break 
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down. In those regions, a refined treatment is deserved. 
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This appendix is devoted to simplify the expressions of B^k, F^, and D^k in case of 
cylindrical symmetric or spherical traps. For a cylindrical trap, the boson and fermion 
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density distribution depend on p 




, z — Rz only. By writing 
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Similarly, one finds the relations: 
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where B^pi F^i and D^p {a, (3 — p or z) take the form, 



m,a;l,(0)iV, (p^), J ^ ^ dp'' dp' 



my^l(0)N>.{^'')J dz dz ' 



Bzp = .2 /r>\ 7vr /-,2\ / pdpdz-^Zp- 



rufcvlfrnf {z^)f J ^ ^ dz dp' 

n 3bf /• , , 977° 2^77° 

-O^z = 2 /n\ Ar / 2\ / Pdpdz—-Z 



mfu;lf{0)Nf {z^)^ J ' ' dz dz 
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For a spherical trap, the density distribution is a function of r = (p^ + z^)^*^^ only. We thus 
define p = rsin6' and z = rcos9, where 9 G [0,7r]. After some straightforward algebra one 
finds, 



F, = F, F, = F, (A4) 



where B, F and D are given by Eqs. (]T6|). 
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Figures Captions 

Fig. 1. The frequencies (upper part) and mixing angles (lower part) of the monopole 
mode in the collisionless regime as a function of the boson-fermion interaction, for a spher- 
ical trap with (a) Qbb = 0.1, (b) 0.5 and (c) 2.11. The other parameters are Nb = Nf = 10^, 
TUb = TUf = m {or Ub = ujf = uo)- The low-lying and high- lying frequency mode are denoted 
by the thick solid and dashed lines, respectively. The analytic solutions obtained by com- 
bining Eqs. (0), ([T8|) and (p2D are also plotted in figures (a) and (b) by thin lines. Note 
that in figure (c), the result is not meaningful at k > k,, ~ 0.69, where the phase separation 
will occur. 

Fig. 2. The same as in fig. 1, but for the quadrupole mode. 

Fig. 3. (a) The boson (solid line) and fermion (dashed line) density distribution for a 
spherical trap with Qbb = 0.5 and k = 1. In this case, the fermions experience a constant 
potential in the region occupied by bosons and therefore uniformly distributed there, (b) 
b = B / K, f = F/k, and d = D/k, as a function the boson-fermion interaction for a spherical 
trap with Qbt = 0.5. Note that precisely at n = 1, B = D = 0. 

Fig. 4. The same as in fig.l, but for the quadrupole mode in the hydrodynamic regime. 

Fig. 5. The frequencies (upper part) and mixing angles (lower part) of the quadrupole 
mode in the hydrodynamic regime against the boson-fermion interaction, for a spherical trap 
with (a) mf/rrib = 0.8 and (b) rrif/mb = 1.03 at Qbb = 0.5. Since the boson and fermion 
trapping frequency is not the same, the mode frequency for bosons and fermions is no longer 
degenerate at k = 0, i.e., it is out of the resonance. As a result, the degree of mixing between 
bosonic and fermion collective oscillations is much reduced. 
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Fig. 6. The frequencies of each mode as a function of k for a cigar-shaped trap with 
A = 0.5. (a) quasi-monopole and (b) quasi- quadrupole. The results in the coUisionless and 
hydrodynamic regime, respectively, are displayed by the thick and thin lines. The other 
parameters are Ni, — Nf = 10^, nib — = m and gi,b — 0.5. 

Fig. 7. The same as in fig. 6, but for a disk-shaped trap with A = 2.0. 
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